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li^; ABSTRACT 

I Context. In the inner parts of the Galaxy the Infrared Dark Clouds (IRDCs) are presently believed to be the progenitors of massive stars and 
^ star clusters. Many of them are predominantly devoid of active star formation and for now they represent the earliest observed stages of massive 
star formation. Their Outer Galaxy counterparts, if present, are not easily identified because of a low or absent mid-IR background. 
Aims. We characterize the ambient conditions by determining physical parameters in the Outer Galaxy IRDC candidate Gl 1 1.80+0.58, a 
relatively quiescent molecular core complex in the vicinity of active star forming regions such as NGC 7538 and S159. 

Methods. We conduct molecular line observations on a number of dense cores in Gil 1.80+0.58. We analyze the data in terms of excitation 
temperature, column and volume density, mass and stability. 

Results. The temperatures we find (15-20K) are higher than expected from only cosmic ray heating, but are comparable to those found in 
massive cores, such as IRDCs. Star forming activity could be present in some cores, as indicated by the presence of warm gas (NH3, '^^CO 
self-absorption) and Young Stellar Object candidates. The observed super-thermal line-widths are typical for star forming regions. The velocity 
^ , dispersion is consistent with a turbulent energy cascade over the observed size scales of the complex. We do not find a correlation between 
J. the gas temperature and the line-width. The LTE masses we derive are much larger than the thermal Jeans mass. Therefore, fragmentation is 
expected and may have occurred already, in which case the observed lines represent the combined emission of multiple unresolved components. 
Conclusions. We conclude that Gil 1.80+0.58 is a molecular core complex with bulk properties very similar to IRDCs in an early, but not 
- - pristine, star forming state. The individual cores are close to virial equilibrium and some contain sufficient material to form massive stars and 
star clusters. The ambient conditions suggest that turbulence is involved in supporting the cores against gravitational collapse, at least down to 
the observed sizes. Additional high resolution data are necessary to resolve and analyze the smaller scale properties. 

Key words. Molecular data ~ Stars: formation - ISM: clouds, molecules, structure 
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1. Introduction 

Infrared Dark Clouds (IRDCs) were discovered a decade 
ago as dark silhouettes against a bright mid-infr ared back- 



grou nd by the Midcourse Space Experiment (MSX, Ega n et al 
1998j) and the Infrared Space Observatory (ISO, Perault et al 
19961). Manv IRDCs contain corapact (sub-) millimeter 



cores jCarev et al."2000t lOarav et al.ll2004l; lOrmel et al.ll2005 



Rathborne et al.M2005i) and the current picture is that these mas- 
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sive dense cores represent the early, cold stages of clus tered star 
formation dRathborne et al.ll2006l:TMenten et al.ll2005i and ref- 
erences therein) and that they are potentially the biith sites of 
massive stars. In order to assess the role of IRDCs in the pro- 
cess of star formation, their ambient physical conditions need 
to be determined. Only then insight may be gained into the dif- 
ferences between low- and high-mass star formation, the nature 
of the initial mass function and the impact of environment on 
star (cluster) formation. That is, the putative early stage that 
IRDCs represent implies that their physical state provides a di- 
rect record of the initial conditions pertinent to stellar birth. 
A complete picture of the star forming properties in massive 
dark clouds, specifically the effects of external conditions on 
the formation process, requires a study of similar objects in 
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different environments. The Outer Galaxy obviously represents 
the place in our Galaxy where the conditions, such as metallic- 
ity, density, interstellar radiation field and overall star format- 
ing activity a re different compared to th e inner regions of the 
Galaxy (e.g.. lBrand & Wouterloolill995l:lRudolph et al.ll2006l) . 
Variations of these properties from cloud to cloud may affect 
the local star formation rate, the mass distribution (IMF) or the 
star forming efficiency. In fact, only a direct comparison be- 
tween Inner and Outer Galaxy star forming regions can asses 
in what way star formation differs with Galactic radius, if at 
all. However, due to a lack of bright background emission else- 
where, IRDCs were identified only in the inner parts of the 
Galax y, mainly toward th e Molecular Ring and the inner spiral 
arms (ISimon et al ]l2006l) . Can we identify massive dark clouds 
in an early evolutionary stage in the Outer Galaxy using a dif- 
ferent approach? 



Frieswijk & ShipmanI (120071 in prep., FS07 hereafter) propose 
a list of candidate IRDC-like objects in the Outer Galactic 
Plane in their investigation of the distribution of highly red- 
dened point sour ces observed in the Two Micron All Sky 
Survey (2MASS, iSkrutskie et al.ll2006l) . They suggest that a 
number of these clustered red sources are reddened due to 
foreground extinction in the form of dark clouds. As a ver- 
ification of the presence of molecular material, they iden- 
tified CO structures, observed in the Five College Radio 
Astrono my Observatory C O Survey of the Outer Galaxy 
(FCRAO. lHever et al.l 19981) . that match the global morphology 
of the regions. Even though CO traces mainly the outer layers 
of molecular clouds, the amount of extinction they find indi- 
cates the presence of large column densities. A subset of the 
structures in this near-IR study are not associated with MSX or 
IRAS point sources and this would also suggest an early star 
forming nature. Additional data, e.g., of molecular lines and 
dust continuum are required to confirm if some of these near- 
IR structures are indeed the Outer Galaxy counterparts to the 
Inner Galaxy IRDCs. 

In this paper we present follow-up molecular line observa- 
tions of the Outer Galaxy dark cloud candidate Gl 11. 80-1-0.58 
in the direction of Cepheus in the Perseus spiral arm. We 
derive physical properties for a number of dense cores lo- 
cated along the filamentary cloud structure, such as tempera- 
ture, density and mass and we characterize the nature of the 
cloud by comparing the results to existing studies. The cloud 
is part of the C as OB2 complex at a radial velocity around 
-55kms-i (e.g..l Blitz et al.ll982l) . There are a few well-known 
star forming regions located nearby, such as NGC 7538 and 
S159 . NGC 7538 is actively forming stars and extensive stud- 
ies have been conducted and reveal various stages of evolu- 
tion including compact dense cores, HII regions, several mas- 
sive proto-stellar obje cts and many lower mass Young Stellar 
Objects (YSO'sl (e.g..lWerner et al.lll979l:lKameva et alJll990t 



Sandell & SieverslbOoiT SI 59 is a bright reflec tion nebulae il- 
lumi nated by an optically visible O-type star dChini & Wink 
[l984) . Part of the cloud is associated with IRAS 23133H-6050 
and seen as an emission nebula. In the immediate surround- 
ing of S159 several compact radio continuum sources and 



Section |2] describes how the object was identified and gives a 
brief overview of the the appearance of the region and its im- 
mediate environment in MSX, 2MASS and FCRAO data. The 
observational setup and the target positions are given in Section 
[3]and Section|4]describes the observational properties. The de- 
rived physical properties of the cores are given in Section |5] 
In Section |6] we discuss the results and present our conclusion 
on the nature of dark cloud Gill .80-1-0.58. Section|7]ends with 
some concluding remarks. 

2. G111 .80+0.58: Source selection 

Gl 1 1.80-1-0.58 was selected from a list of candidate dark 
clouds. This list comprises all sources that were identified in 
the 2MASS Point Source Catalog (PSC) as extended red fea- 
tures in the Outer Galactic Plane (FS07). In this work, a sta- 
tistical measure, using the Mann-Whitney-U-test, is performed 
on a 60" X 60" grid covering the entire Outer Galactic Plane. 
Adjacent cells where the {H - Ks) color distribution of stars is 
different from the local surroundings on a high (> 99%) confi- 
dence level are selected as initial targets. The dark cloud can- 
didates are then chosen on the basis of their abnormal red color 
and an absence of counterparts in SIMBAD. In this Section we 
briefly describe the appearance of the Gl 1 1.80+0.58 region in 
MSX, 2MASS and FCRAO data. Details of interest for the in- 
dividual target positions which were selected for single point- 
ing observations (see Sec. 14.3b are discussed in Section l63] 
We selected this region specifically for a follow-up study be- 
cause it was identified as an excessively red, filamentary region 
in a pilot-survey covering blindly about ten square degrees in 
the second Galactic quadrant. The color scaling in Figures 1, 2 
and 3 shows the MSX 8 jum emission toward the location of the 
complex and its immediate surroundings. The two prominent 
emission regions labeled A and B correspond to NGC 7538 
and S159, respectively. The weaker emission near label C is 
associated with several masers and radio sources. The outlined 
region depicts the area mapped in C'^O (Fig. 4). The target po- 
sitions for single pointing observations (Table 1) are indicated 
by the black open stars. Compared to the active regions in the 
field, only faint 8 fim emission (< 5 x 10"^ Wm^^sr"') is asso- 
ciated with this area. Within the box, several IRAS sources are 
present, indicated by the filled red stars. 
The blue contours in Figure 1 display the red tail of the 
(H - Ks) colors of stars observed in 2MASS on an oversam- 
pled 30"x 30" grid. The solid and dashed contour represent 
the values 0.52 and 0.68 mag, respectively. Assuming that the 
(H - K^) colors are related to the extinction they can be an- 
alyzed using the Near Infrared Color Excess technique (e.g., 
Lada et alJll994l) . The near-IR extinction in a cell is then given 
by 



< E{H - /:s)ceii 



N 



(1) 



(comp act) HII regions have been identified (e.g.. lLebr6n et al 
20011) . 



where the summation is over stars present in a grid-cell. The 
color excess E{H - Ks)i per star is the difference between the 
observed color and the intrinsic color, where the latter is de- 
rived from an off-position chosen to represent the color distri- 
bution free of extinction. Note that foreground stars will play 
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a significant role and reduce the average reddening per cell 
at large cloud distances. A correction for this, e.g., by com- 
paring with the off-position, is required to determine the ex- 
tinction accurately. This is the main reason why the identi- 
fication is based on a statistical color distribution instead of 
the color excess, because a priori a distance to the objects 
is unknown and off-positions were not chosen while process- 
ing the entire Outer Galactic Plane. Keeping this in mind and 
adopting an intrinsic {H - Ks) value of 0.14 mag, the con- 
tour values correspond to ~ 6 and ~ 8.5 mag extinction, for 
convenience converted to Ay using a standard extinction law 
(i.e., Av=15.9x <£■(// - K<,)..n > : l iTeke & Lebofskylll985h . 
The peak extinction that is measured this way corresponds to 
~ 15 mag in Ay. Note that these values are similar to those re- 
ported for Inner Galaxy IRDCs using the same color excess 
method (e.g., G 48 toward W51 with peak Ay ~10-20mag; 
ShipmanetatlboOSl) . 

The contours in Figure 2 give a measure of the 2MASS stellar 
distribution. The mean distribution of the 2 by 2 square de- 
gree suiTounding field is ~ 8 stars per 60" x 60" cell with a Icr 
spread of 3 stars. The red contours depict the regions deficit in 
stars (4 and 3 stars per cell) whereas the blue contours show 
a surplus (11 and 13 stars per cell). The star counts are also 
related to the extinction along the line of sight, but here, a di- 
rect translation is difficult because they represent data from an 
incomplete 2MASS catalog, i.e., the Point Source Catalog in- 
cluding the faint extension. Therefore, they merely give an in- 
dication of the stellar distribution and the number of stars that 
is used for the identification and color excess method. 
Based on the 2MASS data, four different appearances can be 
distinguished in the area of interest and may be explained intu- 
itively by the following: 

1) The combination of color excess and a deficit in star counts. 
The measured extinction is due to background and, if present, 
embedded stars and the column of foreground material is 
large resulting in fewer background stars compared to the off- 
position. 

2) The combination of color excess and a surplus in star counts. 
The measured extinction is due to background and embedded 
sources. The embedded objects result in an enhancement in the 
star count distribution. 

3) Only star count contours. A large column of material pre- 
vents the observation of sufficient background reddened stars 
resulting in a decline of the stellar distribution. The average 
color distribution is normal either due to a complete absence of 
red background stars or due to the presence of too many fore- 
ground stars. 

4) Only color excess contours. Red background and embedded 
sources account for the average red color and the embedded 
sources may compensate for the deficit of background stars due 
to the extinction. 

Note that S159 and NGC 7538 both show up as reddened re- 
gions and have a surplus of stars, as can be expected for these 
active star forming regions where multiple embedded sources 
are present. 

Gl 1 1.80-1-0.58 is not an entirely unknown object. It was iden- 
tified in the FCRAO survey a nd cataloged as part of a molec- 
ular cloud dBrunt et al.ll2003h . The blue contours in Figure 3 




111.8 

1 (degrees) 

Fig. 1: 8 i-im emission of the region toward Gl 1 1.80H-0.58. The 
two prominent sources labeled A and B are NGC 7538 and 
S159, respectively. Region C is associated with several masers 
and radio sources. The region of interest for the current work is 
outlined by the box and is mostly devoid of bright 8 fim emis- 
sion. The contours show the (H - Kg,) color distribution. The 
solid and dashed contours coiTespond to 0.52 and 0.68 mag, 
respectively. The black open stars represent the target positions 
(Table 1) where we conducted single pointing observations. 
The red filled stars correspond to IRAS point sources. 

[10"^ ¥m"^sr"'] 




112 



111,8 111.6 
1 (degrees) 



111.4 



Fig. 2: Same as Figure 1, except the contours represent the 
2MASS source distribution on a 60"x60" grid. The red 
dashed and solid contours correspond to 3 and 4 stars per cell, 
respectively. The blue dashed and solid contours correspond to 
1 3 and 1 1 stars per cell, respectively. The average source count 
per cell is 8. 
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display the CO data integrated approximately between -45 
and -60km s '. The cloud is in velocity space as well as in 




112 111.8 111,6 111.4 

1 (degrees) 



Fig. 3: Same as Figure 1, except the contours represent the 
integrated CO emission between -45 and -60 kms"' as ob- 
served in the FCRAO survey. The contours range from 60 to 
lOOKkms ' in steps of lOKkms '. 

projection on the sky part of a large star forming molecular 
cloud complex in the Perseus spiral arm which also includes 
the aforementioned regions NGC 7538 and S159. The kine- 
matic distance to the Gill. 80-1-0. 58 complex, assuming a flat 
rotation curve and the lAU standard constants V© - 220km s"' 
and/?Q = 8.5 kpc, is ~ 5.0 kpc. S159 has a similar central veloc- 
ity (-56 km s"') as the object presented in this paper However, 
a distance of 3. 1 kpc was reported by Brand & Blitz ( 1993) for 
this cloud, based on spectrophotometric observations of asso- 
ciated stars. Also, the distance to NGC 75 38 is generally con- 
sidered to be less (2.8 kpc; e.g.. Ivan der Tak et a l. 2000) than 
the kinematic distance. It appears that this part of the Perseus 
spiral arm is actually closer to us than suggested by the radial 
velocity. For the analysis in this paper we adopt a distance of 
3.1±0.2kpc, but this may be scaled accordingly. 
More extensive details on the identification of candidate IRDCs 
in the Outer Galactic Plane are given in FS07, but the spectro- 
scopic data presented in this paper indeed confirm the presence 
of a filamentary molecular cloud complex at the identified po- 
sition. 

3. Observations 

The observations presented in this paper were carried out with 
the Effelsberg 100-m telescope in Germany and the IRAM 30- 
m telescope in Spain early 2005. An overview of the observed 
positions and the main parameters of the observed lines are 
given in Tables [1] and |2l 



Table 1 : Source list 



Name 




b 


^1 




n 


n 


(^pc) 


Gill 80+0 58 


111.80 


0.575 




PI 


111.65 


0.59 


1.04 


P2 


111.67 


0.58 


0.80 


P3 


111.67 


0.61 


1.14 


P4 


111.71 


0.59 


0.77 


P5 


111.75 


0.59 


1.01 


P6 


111.76 


0.65 


1.01 


P7 


111.78 


0.50 


0.90 


PS 


111.78 


0.58 


1.19 


P9 


111.82 


0.56 


0.60 


PIO 


111.83 


0.58 


1.13 


S159 (reference) 


111.61 


0.38 


0.60 



" Size of the corresponding core (see Sec. 15. 2t . derived from the 
C'^O map displayed in Fig. 4 and assuming a dista nce £) of 3.1 kpc, 
where D is adopted from S 159 terand & Bhtjl993h . 



Table 2: Line parameters 



Line 




Frequency 


HPBW 


^sys 










(GHz) 


(") 


(K) 






'3CO 


1-0 


110.201 


22.5 


100-200 


0.95 


0.75 




2-1 


220.399 


11.2 


400-900 


0.91 


0.77 




1-0 


109.782 


22.5 


110-200 


0.95 


0.75 




2-1 


219.560 


11.2 


200-600;. 


0.91 


0.77 


c^-^s 


2-1 


96.413 


25.5 


115-135" 


0.95 


0.55 


NH3 


1,1 


23.694 


40 


120-400 




0.58 




2,2 


23.723 


40 


120-400 




0.58 




3,3 


23.870 


40 


120-400 




0.58 



" Forward efficiency, only for IRAM observations 
* Beam efficiency 
200 - 350 K for the map, 400 - 600 K for single pointing observa- 
tions 



3.1. Effelsberg 100-m observations 

Single pointing observations of the NH3 (1,1), (2,2) and (3,3) 
inversion transitions were performed using the Effelsberg 100- 
m telescope of the Max-PIanck-Institut fiir Radioastronomie 
to assess the gas kinetic temperatures at the positions listed in 
Table[T] The first set of data were taken under reasonable winter 
conditions in January and February 2005 using the 1 8 - 26 GHz 
frontend in position switch mode. The receiver was tuned to a 
central frequency of 22 GHz and the auto-correlator with 2 sub- 
units of 20 MHz bandwidth covered the 3 transitions simulta- 
neously. After smoothing the data the spectral resolution was 
K0.2kms At the 100-m telescope, no chopper wheel cali- 
bration is available, so that the atmospherical opacity and the 
gain response of the receiver need to be calibrated with photo- 
metric references. The details of the procedure used to calibrate 
our data are given in Appendix A. The resulting system noise 
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temperatures are of order 120-400K. A beam efficiency of 
- 0.58 is used to convert the antenna temperature to the 
main beam temperature Tmb. Additional data were taken in the 
following months and calibration was performed using a posi- 
tion observed in the first run. Data reduc tion was done using 
the CLASS package (IForveille et al ] ll989h . For the analysis of 
the kinetic temperature the absolute scaling is not a concern be- 
cause it depends only on the ratio of the hyperfine (1,1) lines 
and the (1,1) to (2,2) brightness temperature ratios. The spectra 
are displayed in Figures lBniB.2l and lB.3l in the Appendix. 



3.2. IRAM 30-m observations 

Observations at the IRAM-30m telescope were conducted un- 
der excellent winter conditions (zenith opacity at 220 GHz 
< 0.05) in February 2005. A fully-sampled map of 18'x 18' in 
C'^O 2-1, covering the main c omplex, was obtained using the 
HERA multi-beam instrument (^ Schuster et al. 2004). The cor- 
responding SSB receiver and system temperatures were sjIOO- 
250 K and 200- 350 K, respectively. The on-the-fly observ- 
ing mode was used to allow continuous data acquisition as 



the antenna was moving, with a scanning velocity of l"s" . 
Frequency-switching mode was used to subtract the sky back- 
ground contribution. The HERA matrix projected on the sky 
was rotated by 9.6° to provide a 4" spatial sampling in both 
directions. Maps were done in lambda and beta directions to 
minimize striping eff'ects due to temporal drifts. Each spectrum 
has 896 channels of 80 kHz or 0.1 kms ' width thanks to the 
VESR4 autocorrelator facility backend. Data were reduced us- 
ing the GILDAS softwar^H The bandpass of the system was re- 
moved from each spectrum prior to folding, by subtracting a 
low-order (< 3) polynomial. The spectra were then resampled 
on a 6" grid by Gaussian convolution with a final rms of 0. 17 K 
in each 0. 1 km s ' channel. 

In addition, pointing observations were performed toward the 
integrated intensity peaks in C'^O 1-0 and 2-1 (for column den- 
sity), '^COl-0 and 2-1 (for temperature) and C^'^S2-1 (as a 
high density tracer). The receivers were connected to an au- 
tocorrelator with a resolution of 80 kHz resulting in velocity 
channels of 0.2-0.25 and 0.1 kms"' at 3 and 1.3mm, respec- 
tively. The rms noise levels are of the order of 0.3 and 0.5 K at 
3 and 1.3 mm, respectively. The absolute calibration to obtain 
the main beam temperature Tmb involved a standard conversion 
of the antenna temperature T'^ using the forward and beam ef- 
ficiencies listed in Table|2l The observed intensities, and hence 
all the derived parameters, are overestimated due to a contribu- 
tion of the emission present in the error beam. We adopt a very 
conservative estimate for this contribution of 20%. 



4. Observational properties 

4.1. Spatial distribution: C^^O 2-1 map 

The high-resolution (~ 12") intensity map of the C'^O 2-1 tran- 
sition (integrated between -48 and -56 kms"') is displayed in 
Figure 4. The C'^O molecule, being optically thin under typical 



ISM conditions, traces the gas column along the line of sight 
and reveals the filamentary structure in great detail. Several 
high column density regions are identified along the filaments 
and a number of these were selected for single pointing obser- 
vations (PI -PIO). The extend of the C'**0 emission in the map 
is used to determine the size of the cores (Sec. 15.21 ). 
The reddening contours derived from 2MASS in Figure 1 show 
a similar structure, but they do not trace the C"^0 everywhere. 
This can be explained by both a lower resolution and the limi- 
tations of the color excess method, which relies completely on 
the presence of background and embedded red objects along 
the line of sight. If the extinction is too high, no red background 
sources are detected and the method fails. However, in such a 
case, a decrease in star counts (Figure 2) can be expected and 
indeed this is seen for the regions near PI, P2, P3 and P4. 
Filamentary struct u re is frequently observ ed in IRDCs (e.g.. 



Carey et al.l Il998t iJohnstone et al.l 120031) and is predicted 



by theories of clo ud evolution and clump formation (e.g., 
iKlessen et all2004 . 

In this paper we focus on the physical properties of a num- 
ber of C"*0 intensity peaks, indicated by the open stars and 
corresponding to the positions listed in Table [T] The filled red 
stars represent the location of IRAS point sources in the region. 
Except for P5 and possibly P9 and PIO (see Section |631 ). the 
IRAS sources do not coincide or appear to interact with any of 
the target intensity peaks and though they may be associated 
with the complex, their connection is not considered further in 
the work presented here. 



4.2. Velocity distribution: 2-1 channel maps 



Though it is beyond the scope of this paper to fully analyze 
the velocity maps of the region it is worth showing the spectral 
velocity channels to reveal the complexity and different 
components present in the region. The velocity maps are given 
in Figure 5 and below is a list of characteristic features that 
can be identified: 

fl) the horizontal structure just below the center (e.g., at 
-5 1.5 kms"') containing several of the observed positions 
(P4, P5, P8, PIO); 

b) the curved, vertical structure (e.g., at -53.5 kms"') contain- 
ing position P6 and moving up at nearer velocities; 

c) the main intensity peak (corresponding to P8), possibly 
enhanced in intensity because of the superposition of the 
filamentary structures mentioned above, though the core itself 
is seen already at -51 kms"'; 

d) the positions of PI, P2 and P3 on the right side are dominant 
at the nearest velocities (-48 to -51 kms"') and appear to be 
the end of the horizontal structure, bending towards us; 

e) the structures in the upper left corner, where no target 
positions are located, may be connected with the vertical 
curved structure. 



URL: |http : //www. iraiti. fr/IRMFR/GILDAS/| 
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T„i, [K km s"*] 




111.9 111.8 111.7 111.6 

1 (degrees) 



Fig. 4: This Figure shows the C O 2-1 map of the filamentary dark cloud complex Gl 11.80+0.58 observed with the HERA 
instrument on IRAM, integrated between -56 and -46kms The cores listed in Table 1 are given and correspond to the 
positions indicated by the open black stars, where additional single pointing observations were conducted. The filled red stars 
represent the IRAS point sources. The linear scale is indicated assuming a distance of 3.1kpc. The approximate half-power- 
beam-widths (HPBW) represent the '^CO and C'**0 2-1, the C^^'*S 2-1, '^CO and C'**0 1-0 and the NH3 beams, respectively. 



4.3. Single pointings 

The observed spectra of '^CO, C'^O and C^"*S toward the main 
C'^O 2-1 intensity peaks in the map are given in Appendix iBl 
Furthermore, the Gaussian components that were fitted using 
the GILDAS software are displayed as dashed lines. Some 
positions have additional components in the CO lines that do 
not correspond to a C^'^S detection. The focus of this paper 
however, is on dense cores, i.e., cores that do have a C^"^S 
detection. Therefore, the apparently lower density components, 
likely associated with lower density cores, cloud envelopes or 
specific phenomena such as outflows are not included in the 
analysis. 



We used the following approach to define the dense com- 
ponents; 1) Gauss-fitting the C^"*S with components, 2) 
Gauss-fitting the C'^O and '^CO species with at least lines, 
3) imposing as constraints the positions that correspond to the 
C^'^S fits. The different components are numbered according to 
increasing central velocity (Column 2, Table 3). Some of the 
'^CO lines show high optical depth effects and are impossible 
to match in velocity with the optically thin lines. In that case 
fits are produced with the central velocity as a free parameter. 
Hence the C''^0 is used to determine the column densities and 
masses. 
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Fig. 5: This Figure displays the velocity channels centered on the velocities given in the upper right comers in the range between 
-48 and -55.5 km s"' and integrated over 0.5 km s '. The region shows a complex velocity distribution and several main features 
can be identified (see text). 



5. Results: Physical properties of the clouds 

The physical characteristics of the cores that are determined 
from the molecular lines are listed in Tables 3 and 4. The 
integrated properties (column density and derived parameters) 
as well as the '-^CO temperature are determined by using the 
Gaussian fits of the individual components. A brief description 
of the tables and their content is given here. More details on 
the derivation of the values are presented in the following 
sub-sections: 

Table 3- a) Physical properties assuming the kinetic tem- 
perature given by '■'CO. 

Column (1): Identification number of the cores sorted accord- 
ing to Galactic coordinates. 



Column (2): Identification number n of the component along 
the line of sight; coiTesponds to the n''' component in the 
spectra counting from the left, where n is the same for the 
different molecules if they have a similar central velocity. Note 
that only those components that correspond to C^'^S emission 
are part of the analysis and only these are listed. 
Column (3): Kinetic temperature T]^n (K) derived from '^CO. 
Column (4-5): Column density (in units of 10''*cm"^) of the 
CO molecules assuming LTE and optically thin emission. In 
brackets the optical depth at the line center is given. Note 
that the '^CO emission is in general optically thick in dense 
cores and for some positions only a lower limit could be 
determined for both Tc and Nmoi- 

Column (6): Column density (in units of 10^' cm"^) of H2 
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derived from the C^^O emission u sing a standard abundance of 
1.67 X 10-^ (fFrerkingetal.lll982l) . 



Column (7): Average H2 volume density (in units of 10 cm ) 
derived from the H2 column and using the size. 
Column (8): LTE Mass (M©) given by Equation |9] assuming a 
constant volume density and using the size from Table [T] 

Table 3-b) Physical properties assuming the kinetic tem- 
perature given by NH3. 
Column (1-2): Same as part a. 

Column (3): Kinetic temperature Tkin (K) derived from NH3. 
Column (4-8): Same as part a. 

Table 4) Line-widths and mass characteristics derived 

from the IRAM data. 

Column (1-2): Same as Table 3. 

Column (3-5): Line-width (FWHM) of the observed molecules. 
The '^CO and C''^0 values are an error-weighted mean of the 
1-0 and 2-1 transitions. 

Column (6): The error weighted mean FWHM of the optical 
thin C"*0 and C^"^S lines that is used to derive the virial mass. 
Column (7): Virial mass (M©) as given by EquationfTOlassum- 
ing a constant volume density n(H2) and using the line-width 
Ay given in Column 6. 

Column (8): LTE mass (M©). Same as Column (8), Table 3-a. 
Listed here for comparison with the other masses. 
Column (9): Thermal Jeans mass as given by EquationfTTI 
Column (10): Thermal over gravitational energy ratio 77 as 
given by EquationfT2l 

5.1. Kinetic temperatures 

The kinetic temperature of the cores can be estimated from both 
the NH3 and the '^CO lines. In general, considering the condi- 
tions in dense cores, collisions dominate the excitation process 
and the levels are thermalized. In this case, the kinetic temper- 
ature equals the excitation temperature, i.e., = Tex. A brief 
description of the determination of Tgx for both molecules is 
given below. 

A comparison between the temperatures derived from NH3 and 
'^CO can be misleading because they likely represent different 
regions within the cores. NH3 is a good temperature tracer and 
because the molecules do not fr eeze out onto the dust grains 
for densities below ~ 10'' cm"^ dBergin & Lan ger 1997) they 
are observed even in very dense regions. In the inner, dense 
parts of cores the heating is due to cosmic rays and, if present, 
internal sources. The '^CO molecule is likely depleted in dense 
and cold cores and it is optically thick due to its much larger 
abundance. Thus it traces mostly the outer parts of molecular 
clouds where heating is due to the interstellar radiation field 
penetrating the outer layers. 

5.1 .1 . Ammonia 

The rotational temperatures of the NH3 molecule, characteriz- 
ing the level populations, were determined by fitting the main 
and hyperfine components of the (1,1) transition and the main 



component of the (2,2). The reduction of the (1,1) transition 
was done using "METHOD NH3(1,1)" in CLASS to fit the hy- 
perfine structure, whereas for the (2,2) line a standard Gauss 
fitting procedure was suffic ient. A detailed descr iption of the 
standard analysis is given bv lBachiller et al. ( 1987 ). The kinetic 
temperature is derived using the analytical expression given by 
IXafalla et aL (.2004.) : 



l--^ln[l + l.lexp(-16/r„t)] 



(K), 



(2) 



where Trot is the rotational temperature. Table 5 lists the prop- 
erties derived from the NH3 observations. Columns 1 to 4 give 
the target position, line component, central velocity and line- 
width (FWHM) of the (1,1) transition, respectively. Columns 
5-6 and 9-10 give the main beam temperatures and the noise 
of the spectra for the (1,1) and (2,2) transition, respectively. 
The optical depth r and the corresponding error derived from 
the hyperfine structure fitting for the (1,1) transition are listed 
in Columns 7 and 8, respectively. Typical values for the opti- 
cal depth are around unity with the exception of PI and S159 
CP2. Note however that the error is of the order of r itself and 
the hyperfine fitting is considered dubious. The emission of the 
hyperfine structure for P2, P6 and P9 CP2, if present at all, is 
too weak to be detected at the observed signal to noise lev- 
els. Hence the missing entries. The rotational temperatures are 
given in Column 1 1 with the corresponding error in Column 12. 
Column 13 and 14 give the kinetic temperatures derived using 
Equation|2]and the corresponding error, respectively. Typically, 
the kinetic temperatures range from 13-30K. Note that P2 is 
detected only in the (1,1) transition and P6 not at all, hence the 
missing entry of a central velocity for the latter. The r^b val- 
ues for P2 and P6 are derived from the noise of the observations 
and the temperature parameters merely give an upper limit. 

5.1.2. Optically thick ^^CO 



From basic radiative transfer (e.g., iDickmanI 119781) . the ob- 
served radiation temperature for a certain molecule toward a 
core is given by the expression 



= [/v(7'ex) - Jv(Tbg)][i - exp(-Tv)] , 



in which 



Jv(T) = 



hv 



k exp(^) - 1 



(3) 



(4) 



The constants h, k, Tex, Tbg and Ty are the Planck constant, 
the Boltzmann constant, the excitation temperature, the back- 
ground temperature (~ 2.7 K) and the optical depth at fre- 
quency V, respectively. In the dense cores considered in this 
paper, the '^CO lines are usually optically thick and the results 
derived in Section [53] confirm indeed that Ty exceeds unity at 
the line centers. The radiation temperature of the molecule is 
approximately given by 



r*('^co) = [jy(T,,) - A(rbg)] , 



(5) 



where r^('-^CO) is the observed peak intensity of the '^CO line 
(here we use the 1-0 transition). From this equation the excita- 
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Table 3: Physical properties 

ID CP ru„"cqi auTFj n^. m;7e 

'^COd-O) C'«O(l-0) (i^x (ii^x 

(a) (K) (10"' cm--) (10''' cm-^) 10^' cm-^) 10^ cm-^) Mq) 



(i) 


(2) 


(3) 


(4) 




(j) 




(6) 


(') 


(8) 


S159 


2 


29.1 (±5.1) 


1340 (±751) 


[1.27] 


97 (±31) 


[0.09] 


58 (±19) 


313 (±100) 


204 (±65) 


PI 


2 


10.6 (±1.6) 


>496 


[>4.79] 


11 (±4) 


[0.11] 


6 (±3) 


20 (±8) 


67 (±28) 


PI 


3 


9.8 (±1.4) 


125 (±29) 


[1.44] 


45 (±10) 


[0.41] 


27 (±6) 


85 (±18) 


287 (±62) 


P2 


2 


6.8 (±0.9) 


>195 


[>2.88] 


23 (±5) 


[0.39] 


14 (±3) 


56 (±13) 


87 (±20) 


P2 


3 


10.2 (±1.5) 


589 (±228) 


[4.38] 


24 (±6) 


[0.22] 


15 (±4) 


59 (±16) 


91 (±24) 


P3 


1 


8.7 (±1.3) 


107 (±26) 


[2.02] 


9 (±4) 


[0.14] 


5 (±2) 


15 (±6) 


66 (±29) 


P3 


2 


13.1 (±2.1) 


656 (±547) 


[4.51] 


38 (±9) 


[0.20] 


23 (±5) 


65 (±15) 


291 (±67) 


P4 


1 


7.5 (±1.0) 


51 (±13) 


n 301 


16 (+4) 


[0 271 


10 (+3) 


41 (+11) 

-r J. \^ — ^/ 


57 (±15) 


P4 


2 


10.0 (±1.4) 


>488 


[>3.68] 


21 (+5) 


[0 181 


13 (+3) 


53 (±13) 


73 (±18) 


P5 


2 


9.7 (±1.4) 


227 (±36) 


[1.84] 


26 (±5) 


[0.19] 


15 (±3) 


50 (±10) 


154 (±31) 


P6 


1 


12.4 (±1.9) 


168 (±53) 


[1 541 


17 (+7) 


[0 151 


10 (+4) 


32 (±13) 


99 (+41) 


P6 


2 


12.1 (±1.9) 


>655 


[>4.60] 


48 (±10) 


[0.28] 


29 (±6) 


93 (±20) 


289 (±62) 


P7 


2 


13.2 (±2.1) 


199 (±57) 


[1.24] 


24 (±7) 


[0.14] 


14 (±4) 


52 (±16) 


114 (±34) 


P7 


3 


10.2 (±1.5) 


>318 


[>5.13] 


18 (±8) 


[0.28] 


11 (+5) 


39 (+17) 


86 (±38) 


P8 


2 


13.9 (±2.4) 


>625 


[>3.69] 


124 (±26) 


[0.54] 


74 (±16) 


202 (±43) 


1025 (±217) 


P9 


2 


9.4 (±1.3) 


67 (±18) 


[1.00] 


10 (±3) 


[0.11] 


6 (±2) 


33 (±11) 


22 (±7) 


P9 


3 


13.0 (±1.9) 


250 (±88) 


[2.29] 


11 (±5) 


[0.10] 


6 (±3) 


34 (±17) 


22 (±11) 


PIO 


3 


11.2 (±2.0) 


>389 


[>2.29] 


46 (±11) 


[0.36] 


28 (±7) 


80 (±20) 


347 (±85) 






rki„(NH3) 




















(K) 
















S159 


2 


29.3 (±3.4) 


1341 (±742) 


[1.25] 


97 (±25) 


[0.09] 


58 (±15) 


315 (±82) 


205 (±54) 


PI 


2 


14.4 (±2.6) 


434 (±108) 


[1.50] 


12 (±6) 


[0.07] 


7 (±4) 


23 (±11) 


76 (±38) 


P2 


2 


19.0 (±8.4) 


149 (±39) 


[0.36] 


30 (±20) 


[0.08] 


18 (±12) 


73 (±48) 


113 (±74) 


P5 


2 


18.2 (±2.7) 


217 (±46) 


[0.45] 


34 (±9) 


[0.08] 


20 (±5) 


65 (±17) 


201 (±53) 


P8 


2 


16.7 (±1.4) 


470 (±159) 


[1.63] 


131 (±29) 


[0.40] 


79 (±17) 


215 (±47) 


1090 (±240) 


P9 


2 


20.0 (±2.3) 


78 (±27) 


[0.26] 


15 (±7) 


[0.04] 


9 (±4) 


47 (±24) 


31 (±16) 


P9 


3 


18.0 (±2.6) 


219 (±92) 


[0.90] 


13 (±8) 


[0.06] 


8 (±5) 


41 (±26) 


26 (±17) 


PIO 


3 


13.5 (±2.1) 


435 (±141) 


[2.68] 


49 (±13) 


[0.07] 


29 (±8) 


84 (±22) 


366 (±94) 



" The kinetic temperature is assumed equal to the excitation temperature. 

* Part b is the same as part a, except here the kinetic temperature is derived from NH3. Note that the results are very similar. 



tion temperature can be derived: 

5.29 



ln[l + ™ 



5.29 



(K). 



(6) 



(i3CO)+0.8 



The values are listed in Column (3) of Table 3-a. 
In the following analysis the '^CO 1-0 temperature is adopted 
as the kinetic temperature of the clouds. When the data allow, 
the analysis is also performed using the kinetic temperature de- 
rived from the NH3 lines. 

5.2. Core sizes 

We define a size for each core as the average of the minor 
and major axis of the 50% peak intensity level observed in the 
C'**0 2-1 map. The sizes are given in parsecs in Column 4 of 
Table [1] assuming a distance to the cores of 3.1 kpc. Several 
positions have multiple components along the line of sight and 
obviously the size is then determined from a superposition of 
several cores. However, a distinction between the components 
is difficult because of the small separation in velocity so that 



the same size is adopted for all components. Position P9 is not 
identified as a coherent core in the map and a size of 40" is 
adopted. P4 is part of a filamentary structure and the size is 
assumed equal to the projected width of the filament. 



5.3. Column densities 



The following approach (see e.g., iBourke et al.ll 1997ft assumes 
these conditions: a) the cores are isothermal, and the kinetic 
temperature, Tkin, is given by the excitation temperature de- 
rived from '^CO or NH3; b) the cores have a constant density; 
c) the cores are in local thermodynamic equilibrium (LTE). 
Then, using the observed C'^O radiation temperature (K), 
the optical depth can be calculated as a function of the LSR 
velocity V (km s"') parameterized by the frequency v: 



Tv(C'^0) = -InJl-^ 



1 



kJy(Tbg) 



exp(^)-l 



hv 



(7) 
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Table 4: Line-width and mass properties 

CP Av M~ m;;^ m;; 

'3 CO C'*0 C^'^S virial 







(kms"') 


(km s"') 


(km s"') 


(km s"') 


(tt%;Mo) 




(Mo) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


S159 


2 


2.84 (+0.01) 


2.40 (±0.02) 


2.25 (±0.04) 


2.36 (±0.02) 


352 (±59) 


204 (±65) 


15 (±5) 


0.09 (±0.03) 


PI 


2 


3.00 (+0.06) 


1.07 (±0.15) 


1.57 (±0.25) 


1.22 (±0.13) 


162 (±38) 


67 (±28) 


13 (±4) 


0.17 (±0.08) 


PI 


3 


1.77 (+0.08) 


1.74 (±0.05) 


2.27 (±0.50) 


1.77 (±0.06) 


342 (±40) 


287 (±62) 


5 (±1) 


0.04 (±0.01) 


P2 


2 


2.19 (±0.02) 


1.95 (±0.03) 


1.33 (±0.33) 


1.91 (±0.04) 


307 (±40) 


87 (±20) 


4 (±1) 


0.07 (±0.02) 


P2 


3 


3.98 (±0.02) 


1.62 (±0.02) 


1.67 (±0.85) 


1.62 (±0.03) 


220 (±29) 


91 (±24) 


7 (±2) 


0.09 (±0.03) 


P3 


1 


1 34 C+O 02) 


1.17 (±0.04) 


1 47 (+0 29) 


1 20 C+0 05) 


172 (±20) 


66 (■+29) 


11 (+3) 


16 f+0 07) 


P3 


2 


2.91 (±0.02) 


2.07 (±0.03) 


1.06 (±0.34) 


2.02 (±0.03) 


489 (±46) 


291 (±67) 


10 (±3) 


0.05 (±0.01) 


P4 


1 


1.28 (±0.02) 


1.63 (±0.05) 


0.82 (±0.19) 


1.51 (±0.05) 


185 (±27) 


57 (±15) 


5(±1) 


0.10 (±0.03) 


P4 


2 


2.52 (±0.02) 


1.74 (±0.05) 


1.26 (±0.12) 


1.63 (±0.05) 


215 (±30) 


73 (±18) 


7 (±2) 


0.11 (±0.03) 


P5 


2 


2.83 (±0.02) 


2.25 (±0.04) 


2.62 (±0.24) 


2.29 (±0.05) 


556 (±59) 


154 (±31) 


7 (±2) 


0.07 (±0.02) 


P6 


1 


1.49 (±0.01) 


1.14 (±0.03) 


0.96 (±0.16) 


1.12 (±0.03) 


132 (±15) 


99 (±41) 


13 (±4) 


0.13 (±0.06) 


P6 


2 


2.88 (±0.01) 


2.09 (±0.02) 


1.20 (±0.01) 


2.08 (±0.03) 


459 (±47) 


289 (±62) 


7 (±2) 


0.04 (±0.01) 


P7 


2 


1.91 (±0.02) 


1.49 (±0.11) 


1.62 (±0.30) 


1.52 (±0.11) 


218 (±40) 


114 (±34) 


11 (±3) 


0.11 (±0.04) 


P7 


3 


1.41 (±0.01) 


1.05 (±0.02) 


1.68 (±0.77) 


1.06 (±0.02) 


106 (±13) 


86 (±38) 


9 (±3) 


0.11 (±0.05) 


P8 


2 


2.32 (±0.02) 


2.22 (±0.01) 


2.08 (±0.05) 


2.20 (±0.01) 


604 (±51) 


1025 (±217) 


6 (±2) 


0.02 (±0.00) 


P9 


2 


1.55 (±0.02) 


1.48 (±0.08) 


1.57 (±0.22) 


1.50 (±0.08) 


142 (±28) 


22 (±7) 


8 (±2) 


0.27 (±0.10) 


P9 


3 


1.64 (±0.02) 


1.23 (±0.06) 


1.82 (±0.28) 


1.31 (±0.07) 


109 (±21) 


22 (±11) 


13 (±5) 


0.36 (±0.19) 


PIO 


3 


2.22 (±0.71) 


1.60 (±0.02) 


1.99 (±0.08) 


1.66 (±0.02) 


326 (±30) 


347 (±85) 


7 (±2) 


0.04 (±0.01) 



Table 5: NH3 derived properties 



ID 


CP 


Vlsr 


AV 


r„,b(l, 1) 


Er,„b(i.i) 


T(l,l) 


Er(I.l) 


r„,b(2,2) 


Er„b(2.2) 














(km s-') 


(kms-') 


(K) 


(K) 






(K) 


(K) 


(K) 


(K) 


(K) 


(K) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


S159 


1 


-58.01 


1.13 (±0.17) 


0.59 


0.06 


1.3 


0.7 


0.19 


0.05 


17.1 


2.0 


18.2 


2.0 




2 


-55.87 


2.16 (±0.25) 


0.45 


0.06 


0.1 


0.1 


0.32 


0.05 


25.7 


3.4 


29.3 


3.4 


PI 


2 


-50.12 


2.15 (±0.19) 


0.74 


0.10 


0.3 


0.3 


0.13 


0.07 


13.7 


2.6 


14.4 


2.6 


P2 


2 


-50.95 


1.50 (±0.00) 


0.41 


0.20 






0.14 


0.14 


17.6 


8.4 


19.0 


8.4 


P5 


2 


-52.76 


2.25 (±0.05) 


1.58 


0.09 


1.1 


0.1 


0.50 


0.19 


17.0 


2.7 


18.2 


2.7 


P6 








0.12 


0.12 






0.08 


0.08 


26.4 


23.9 


30.2 


23.9 


P8 


2 


-52.83 


2.07 (±0.08) 


2.11 


0.15 


1.0 


0.2 


0.55 


0.12 


15.8 


1.4 


16.7 


1.4 


P9 


2 


-54.18 


2.82 (±0.82) 


0.54 


0.06 






0.21 


0.05 


18.5 


2.3 


20.0 


2.3 




3 


-52.03 


1.40 (±0.82) 


0.47 


0.06 


1.0 


0.1 


0.14 


0.05 


16.8 


2.6 


18.0 


2.6 


PIO 


3 


-52.23 


1.65 (±0.11) 


1.35 


0.14 


0.8 


0.3 


0.20 


0.10 


13.0 


2.1 


13.5 


2.1 



The corresponding column density can be estimated using 



A^(C'^O) = 2.42 X 10'''^ 



(cm^O , (8) 

l-exp(-^) 

where /i, v, k, JviTh%) are the same as before and AV is the 
velocity step size in kms ^ The column densities derived for 
C'**0 1-0 are given in Column (5) of Table 3. The values listed 
in brackets corresponds to the optical depth, Tc, at the line cen- 
ter and are typically much less than unity for C"^0. 
An analogous derivation is done for the '^CO 1-0 transition and 
these values are listed in Column (4). However, as noted be- 
fore, in the dense cores the '^^CO is usually optically thick (i.e., 
Tc > 1) and the lines cannot be used to determine a reliable 
column density. Nonetheless, the results are listed and for the 
cases where the optical depth becomes too large, a lower limit 
is given. The total molecular column densities, A^(H2), are de- 
rived from the C^^O 1-0 results using a canonical abundance 



of 1.67 X 10"^ ( jFrerkingetaIlll982h . These values are given in 
Column 6 of Table 3. 



5.4. Volume densities 



Column 7 (Table 3) lists the volume density derived from the 
H2 column density, assuming the depth of the core equals the 
size. The density is by definition an average over the actual 
density profile along the line of sight. In the inner parts of the 
cores the density is considerably higher, as suggested by the 
presence of C^'^S. 
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5.5. Mass estimates 

5.5.1 . LTE mass 

Given the size and adopting a constant volume density, the 
molecular mass derived from the lines, referred to as the LTE 
mass, is determined using 

Mlte = A, me «(H2) ( ^5 )3C (Mo) , (9) 

where n{H2) is the volume density, S is the size, niu is the mass 
of a hydrogen atom, fi-2.33 is the mean molecular weight 
consistent with a 25% mass fraction of helium and C = Ipc^ 
lMo~ 1.48x10^^ cm^/g, the conversion to solar masses. The 
values are listed in Column 8 (Table 3). Because the volume 
density is derived from the column density (n(H2) ~ N(li2)/S ), 
the LTE mass of the cores scales with the square of the size and 
thus the distance (i.e., Mlte ~ and ~ D^). 

5.5.2. Virial mass 

Comparing the LTE mass with the virial mass enables an evalu- 
ation of the dynamical state of the cores. Assuming a spherical 
core with uniform density and considering only thermal and 
dynamical broadening (neglecting e.g., magnetic fields, inter- 
nal he ating) the virial mass c an be derived using the expression 
(e.g., iMacLaren et"ai]|l988h 



Mvir = 210/?(Ay)2 (Mo), 



(10) 



where R (=f ) is the radius of the core in pc and AV is the 
FWHM of the line in km s ' . In our case a weighted average of 
the observed optically thin lines, given in Column 6 of Table 4, 
is used. The resulting virial masses are given in Column 7 of 
Table 4. Equation[TO]assumes a constant density. The constant 
may be replaced by 190 or 126 for a density profile given by 



where Mlte is given in Column 8 (Table 4), T is the ki- 
netic temperature listed in Column 3 (Table 3-a), S is the size 
of the core, described in Sec. 15.21 mn, k and G are the 
mean molecular weight (- 2.33), the mass of a hydrogen atom, 
the Boltzmann constant and the gravitational constant, respec- 
tively. Typically, rj is much less than unity and it suggests that 
the thermal energy by itself is insufficient to support the cores 
against gravitational collapse. 

6. Discussion and conclusion 

What kind of object is Gl 1 1.80-^0.58? 

Before we attempt to answer this question it is worthwhile to 
discuss some key properties derived from the observations and 
compare the results with other studies. 



6.1. Gas temperature 

The gas temperatures derived from NH3 (13-20K) and, 
less obvious also from '^CO (7-14K), are higher than ex- 
pected inside molecular cores if only cosmic ray ionization 
is considered as a heating source (a mean Galactic value 
of 4"c R ~ 3 X 10~'^ s~' yields a te mperature of T ~ 8 - 



r ' or p 



r , respectively. 



lOK, Ivan der Tak & van Dishoeckl 12000.) . The NH3 temper- 
atures suggest that deep inside the cores additional heating 
sources are present. Because the '^CO is usually optically thick 
at the observed positions, it traces the temperature of the ma- 
terial at an optical depth of about unity. The actual location of 
this T-l surface in the cores cannot be determined, but it is 
likely that the two molecules trace the temperatures at different 
depths in the cores. 

The derived temperatures are in agreement with values 
found for massiv e and dense cores (e.g., <20K for IRDCs, 
Carev et al.lll998t) . The temperature found for S159 is signif- 
icantly higher (29 K), which is consistent with its more ad- 
vanced star forming state. 



5.5.3. Jeans mass 



6.2. Core masses and total mass 



Thermal Jeans fragmentation, often referred to as the mec ha- 
nism to set the mass scale for star formation ( Larson|[l985 ). is 
approximately given by 



90 



(Mo) 



(11) 



where T is the gas temperature and n is the total volume den- 
sity. The values are listed in Column 9 of Table 4. The typical 
value is of the order of 10 Mo, much smaller than the derived 
LTE masses of the cores. Therefore, structure on much smaller 
scales may be expected within the cores. 

5.5.4. Stability 

The energy ratio listed in Column 10 (Table 4) is defined as 
the thermal energy over the gravitational energy: 



■-grav 



3 ^^LTE 

3 GM?^, 



(12) 



The individual core masses (Mlte) vary from about 20 to 
1000 Mo (approximately IOO-6OOM0 for Mvk) with an av- 
erage of around 190 Mo (about 280 Mq for Myir). Most cores 
are massive when compared to low-mass dark clouds (e.g., 
Evansl[l999 . and references therein), but correspond well with 



masses found for hig h-mass protostellar objects (HMPOs; e.g.. 



Wilhams et al.l 120041) . The LTE mass we find for S159 is in 



agreement with previous high resolution studies ( iLebron et al 



i2001.) where the core is resolved into multiple substructures. 
A rough estimate for the total mass of the complex, includ- 
ing only the observed cores in this paper, is 3000 Mo (about 
4500 Mo for Myh). This i s comparable to m asses found for 
Inner Galaxy IRDCs (e.g., ICarevet al.ll2000l) . 



5 0.5 s 



6.3. Physical characterization of the cores 

The observations suggest that the thermal energy is insufficient 
to support the cores against gravitational collapse. Additional 
support is seen in the observed line-widths toward the cores. 
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These are much broader than expected from purely thermal 
motion (see Column 3-5, Table 4 and Column 4, Table 5). The 
average line-width is about 1.8 + 0.4 kms"'. The sound speed 
AVs, considering that the cores have a typical temperature T of 
1 5 - 20 K, is given by 

/81rKW\^^,,^^_,^ (13) 
and thus the velocities are supersonic. 

The broad lines could be due to turbulent motions in the clouds 
and are indicative of r egions of high-mass star formation (e.g., 
Mvers & Fuller! 119931) . in contrast to mostly thermal broad- 



ened lin es often seen toward low-mass pre-stellar cores (e.g., 
[Ooodma n et al.l 1998 ). The line-width across the filament, ex- 
tending over lOpc (at a distance of 3.1 kpc), is about 5 kms"'. 
The relation between velocity dispersion of substructures (AV/) 
and the velocity dispersion over the larger complex (AVl) in a 
turbulent medium can be expres sed using the empirically de- 
rived Larson's law ( Larson|[l981 ): 



AV, ~ ^Vl 



0.38 



(14) 



where L is the size of the filament and / the size of the sub- 
structure. For a typical core size in our sample of 1 pc, the 
expected line-width is then si2kms ' for individual cores. 
The values we find are consistent with this relationship. The 
standard relations of mass versus line-width are, however, for 
much smaller cores and the larger cores presented here may 
be a superposition of such small cores. Extending the rela- 
tion to smaller cores of 0.1 pc, the expected line-width using 
Equation [T4l is about 0.8 to 0.9 km s"'. This is still in agree- 
ment with values found for typical massive starless cores (e.g., 
Caselli & Mversl[T9 95). 



A decay of the turbulence is exp ected on short timescales (see 



^e 



Ballesteros-Paredes et al. 20061 for a review on star forma 



tion and turbulence) and some mechanism is needed to sustain 
the turbulent support. Internal sources, e.g., deeply embedded 
YSO's may be present and stir up the material from the inside 
through outflows. Indicative of the presence of a heating source 
can be the high temperatures traced by NH3. However, we find 
no correlation between the line-width and the gas temperature. 
The line-width could also be the result of a number of small 
dense clumps at slightly different velocities within the spatial 
resolution element of the observations. This idea is supported 
by the presence of multiple components along the line of sight 
toward many of the observed positions. In addition, the de- 
rived Jeans mass (Column 9, Table 4) for the cores is typically 
~ IOMq and therefore, fragmentation is expected to occur on 
smaller mass scales compared to the observed core masses, if 
thermal processes dominate stability. 

Figure 6 shows a plot of the virial mass versus the LTE mass. 
The virial parameter a can be defined as the ratio of the two 
and a value of unity corresponds to a virial equilibrium state 
(dashed line). The observed values for a range from about 1 to 
3 and thus suggest that, on a global scale, the cores are close 
to virial equilibrium. Note that depletion of the C"^0 molecule 
may play an important role for densities in excess of 10"^ cm"^ 



600 



400 



200 




400 600 
Mlte (Wo) 



Fig. 6: The virial mass is plotted against the LTE mass. The 
dashed line corresponds to equal masses, where the virial pa- 
rameter a - \. The solid line is the fit to the data points. 



if T^dust ~ Tgas (iBergin et al.l 119951) and consequently, the LTE 



mass could be underestimated. Core P8 is an exceptional case 
where Mlte > M^i^. As seen from the channel maps (Fig. 5, 
Sec. 14.2b . in projection P8 appears to be a superposition of two 
crossing filaments. The material traced by the C"*0 may just be 
the line of sight sum of two components. Alternatively, because 
the components have a similar LSR velocity, a collision may be 
taking place between the two filaments. This would lead to the 
formation of a dense (perhaps gravitationally unstable) central 
core. If so, then we may be witnessing this scenario prior to the 
system being fully virialized. 

6.4. Star forming activity 

The NH3 lines are very sensitive to the temperature and there- 
fore the presence of the higher transitions, particularly of the 
(3,3) line, is a good indicator for the occurrence of warmer 
gas, possibly indicating internal heating sources. The (3,3) 
line is clearly seen toward S159, as expected from work by 
Lebron et al.l (1200 lb . and also position P5, where it may result 
from the presence of the associated IRAS point source. A very 
weak detection is seen toward positions P8 and P9. 
Signatures of self-absorption in the '"^CO lines can be an addi- 
tional indication of the presence of the warmer gas inside the 
cores. This may explain also the failure in fitting the '^CO lines 
at the central velocities of the optically thin lines. Most notice- 
able are the positions toward S159, P5 and P8. 
If embedded sources are present, they may have been identified 
in the 2MASS data. Figure 7 shows the (H-Ks) versus (J-H) 
colors of the sources present in the 2x 2 square degree field. 
The grey scaling represents the overall distribution of the field 
and the majority of the sources are located near the dwarf and 
giant sequence, indicated by the black solid curves. The spread 
in the colors is significant but can be addressed entirely by the 
use of the very faint sources in the 2MASS catalog (the faint 
extension) where the photometric errors are larger. The das hed 
lines give the reddening vectors (.Rieke & Lebofskvll 1985b for 
the dwarf and giant sequence and the arrow corresponds to an 
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Ay of 10 magnitudes. In addition, the figure shows for each of 
the cores the associated 2MASS sources in different symbols. 
Note that most star colors toward the cores resemble main se- 
quence stars with or without reddening, and these are likely 
reddened background or 'normal' foreground stars. However, 
the very red (H - Ks) area right of the reddening vectors corre 



spend s to colors of proto-stellar objects (e.g., Robitaille et al 



20061) . Most promising candidates for YSO's correspond to 
sources associated with S159, P5 and P8. 




3.70 



Fig. 7: 2MASS color-color plot. The grey-scales represent the 
location of colors of all stars in the 2 by 2 square degree field. 
Lighter colors correspond to a higher occurrence. The major- 
ity of the stars are situated with some spread around the dwarf 
and giant sequences overlaid as black curves. The arrow cor- 
responds to a visual extinction of 10 magni tudes, in the direc- 
tion o f the reddening vectors (dashed lines, iRieke & Lebofskv 
1985b . The sources associated with the cores are plotted us- 
ing the symbols listed in the legend. The number of associated 
sources is given in brackets. 



6.5. Notes on individual cores 

PI, P2, P3 and PI: These are regions of the molecular cloud 
where multiple components with enhanced densities (based on 
C"^^S) appear to be present along the line of sight. There are no 
clear indications for star forming activity from '^CO self ab- 
sorption, 2MASS colors or NH3 heating (NH3 is only available 
for PI andP2). 

P4: This is part of the horizontal filamentary structure. No clear 
core is seen in the C'**0 2-1 map but high density gas is present, 
indicated by the C^"^S line. Gravity may have induced a density 
enhancement along the filament. Two components are fitted by 
the software to the spectra, but possibly opacity effects cause 
a non-Gaussian behavior of a single component. Manual fit- 



ting of one component would result in an LTE mass of about 
138 Mq, a velocity width of 1.8kms"' and a virial mass of 
262 Mq. 

P5: This position may be associated with, and heated by, the 
nearby source IRAS23137H-6105, which can explain the pres- 
ence of the higher NH3 transitions. One object in the 2MASS 
data has typical colors of a YSO. Also here, the spectra suggest 
that a single manual fit may be favored to two components. 
This would result in an LTE mass of about 220 Mq, a velocity 
width of 2.57 km s"' and a virial mass of 702 M©. 
P6: Two dense (C^"^S) components are present along the line 
of sight. The stellar density appears slightly enhanced but this 
may be a chance encounter along the line of sight. There is no 
indication for the presence of YSO's and none of the NH3 tran- 
sitions is detected which suggests a quiescent, cold stage. 
P8: This is the most profound core in the region and the best 
candidate for being an object similar to Inner Galaxy IRDCs 
(see Section |676] |. Automatic fitting of the CO lines resulted in 
two components, but in C'''*S, a single profile is seen. Manual 
fitting of a single component to the CO lines results in an LTE 
mass of about 1250 M©, a velocity width of 2.28 km s"' and a 
virial mass of 647 M©. This makes the case for being an Outer 
Galaxy IRDC even stronger. Star formation may have started 
here, indicated by the NH3 lines and the presence of YSO can- 
didates in the 2MASS color plot. Opacity effects such as self 
absorption due to colder foreground gas can explain the non- 
Gaussian appearance of the CO lines. The stellar density de- 
rived from 2MASS is significantly higher (1 1-13 stars per cell), 
compared to the rest of the field. 

P9,10: There is arc-like MSX 8/im emission nearby these cores 
(see Figures 1-3). This may be caused by heating due to ra- 
diation from the neighboring source 1RAS23143+6103. Such 
heating may also explain the detection of NH3 (2,2) and faint 
NH3 (3,3) emission. Compression of the material due to a 
shock may then lead to the observed high densities. No indi- 
cation of embedded sources is seen from the 2MASS data. 



6.6. Conclusion: Is G111 .80+0.58 similar to Inner 
Galaxy IRDCs ? 

The global characteristics of the Gil 1.80+0.58 complex that 
are presented in this paper, such as its size (~ lOpc), col- 
umn density (peak> lO^^cm"^) and LTE mass (~3OOOM0) 
indicate a clear resemblance to th e bulk properties of Inner 
Galaxy IRDCs jSimon et alj|2006l) . However, except for P8, 
the individual core properties (e.g., size: ~1 pc, LTE mass: 20- 



350 Mq, density: lO"* - 10"* cm"^) fall short of the values found 
for compact (sub)mm cores in IRDC s (size: 0.02- 0.8 pc, m ass: 



Rathborne et"aDl2006h . 



10- 10-^ Mo, density: lO"*- 10' cm- 
Would this cloud be seen in extinction were it observed toward 
a bright mid-IR background? 

A typical column density found for the cores is 10- 
20xl0^'cm-^, which corresponds to an Ay of 5- 10 mag 



(IBohlin et al.lll978h . The corresponding extinction at 8 /im can 
be expressed as As^,m ~ 0.04Av < 0.4- mag, using conversions 



adopte d from llndebetouw et al.l ( 120051) and iRieke & Lebofskv 



(Il985h . Typical extinctions around Sfim for Inner Galaxy 
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IRDCs are 1-2 mag dCarev et al. I l2000h . Thus, based on this 
simple reasoning the cloud would not meet the IRDC crite- 
rion of mid-lR extinction. However, following the same analy- 
sis, the extinction toward position P8 (74 x 10^' cm"^) is about 
1 .4 mag at 8 fim. Considering that the above values represent 
the densest component of the two CO fits toward PS, the ac- 
tual 8 fim extinction may even be higher when using a single 
fit (1.8 mag). This central part of the cloud makes the region 
a very promising Outer Galaxy IRDC candidate. The presence 
of possible embedded heating sources at this position may in- 
dicate that the core is not in a very ea rly, quiescent state, but 



already forming stars. Recent studies (IRathborne et al.l 12005 
[Ormel et al. 2005, v.d. Wiel, privat comm.), using submm 
and Spitzer data indicate that some Inner Galaxy IRDCs are 
presently forming stars in their cores and along the filaments as 
well. 

Obviously, the dark cloud candidate presented in this paper is 
in the vicinity of star forming activity (NGC 7538, S159) and 
in that sense the boundary conditions may not be that deviant 
from Inner Galaxy IRDCs. It should be noted however, that 
there are environmental difference s in the Outer Galaxy (e.g. , 
radiation field, densi ty, abundance; Brand & Wouterloogll995 



Rudolph et al. mod and references therein) when compared 



to the inner spiral arms and the Molecular Ring, where most of 
the IRDCs are found. 

The super-thermal line-widths and the presence of massive, 
cold and dense cores in Gill. 80-1-0. 58 are both indicative of 
a high-mass star forming complex. The combined mass of the 
cores is perhaps too low to form a massive star cluster (like 
Orion). However, this region is part of a much larger molecular 
cloud complex. We conclude that the Gl 1 1.80H-0.58 complex 
likely belongs to a category of objects simil ar to intermediate- 
mass IRDCs (e.g., IRDC G48 toward W51; lOrmel et al.ll2005l) 
in an early, but not pristine stage of star formation. 

7. Concluding remarks 

Radiative transfer codes can provide a better understanding of 
the properties of the gas, e.g., radial density and temperature 
profiles can be investigated, but require more detailed obser- 
vations. We did co nduct tests with a multi-zone es cape prob- 
ability code (p 3D. IPoelman & Spaansll2006l lioosi) and these 
resulted in similar temperature, density and mass estimates as 
presented in this paper We therefore conclude that for the in- 
tegrated properties these codes have no additional contribution 
to the results. The main advantage when using these models 
is getting a better understanding of the line profiles. In par- 
ticular, effects such as self absorption or velocity structures, 
like infall motion may explain some of the double peaked and 
non-Gaussian features (see e.g., Evans 1999). In the work pre- 
sented here, non-Gaussian features are treated as separate com- 
ponents. Note in this respect that different optically thin lines 
(e.g., C'**0 and C'''*S for core PI) show their peak emission at 
opposite sides of the central velocity Vlsr- 
High spatial and spectral resolution in addition with radiative 
transfer models are required in the future to get better con- 
straints on the physical properties on smaller scales. A spa- 
tial resolution at sub-arcsecond scales may resolve the pres- 



ence of small high density clumps. Resolving this substructure 
will allow an investigation of the Dense Clump Mass Function 
(DCMF). It is of great importance to conduct these studies to 
put constraints on the origin of the stellar mass spectrum. Not 
only in nearby, generally low mass s tar forming clouds (e.g.. 
Pipe Nebula; Lombardi et al. 2006h . but particularly in the 
more massive and distant IRDCs, including the Outer Galaxy 
region presented in this paper 

A different way of characterizing IRDCs is by using (sub)mm 
continuum observations. All In ner G alaxy IRDCs show strong 
submm emission (e.g.. ICarev et al.l 120001) and some contain 
bright, centrally peaked cores. In addition, these data allow 
an independent measurement of properties such as tempera- 
ture, mass and luminosity. It will be worthwhile to observe the 
Gl 1 1.80-1-0.58 cloud in dust continuum and compare the data 
with the molecular line results presented here and with existing 
submm studies in the Inner Galaxy. 
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Appendix A: Calibration of data obtained at the 
100-m telescope 

At the frequency considered here (23 GHz), the Effelsberg tele- 
scope offers a calibration based on noise diodes, so that the re- 
ceived signal is not converted into K but into internal receiver 
counts. In order to both convert the receiver counts into an- 
tenna temperature, as well as to correct for the atmospherical 
absorption, we have observed a photometric calibrator (here 
NGC 7027) at various elevations in the course of each observ- 
ing night. This approach is similar to the so-called "antenna 
tipping" used to measure the sky opacity during e.g., bolome- 
ter observations. The weather conditions were stable enough 
to assume that the opacity derived through this technique was 
reasonably representative of the daily set of data. If Fim is the 
conversion factor between internal counts and Jy above the at- 
mosphere and corrected for backward losses, the measured re- 
ceiver counts are (in a given spectrometer channel): 



r™ = G(el) X e' 



(-TzA) 



, X- 



(A.l) 



where A is the airmass at elevation el, r, the zenith opacity and 
G(el) the normalized elevation gain. Observing a photometric 
calibrator of known flux 5v,ief one has: 
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Using the calibrated flux published bv lOtt et alj (119941) and ob- 
serving the reference at various elevations, the calibration pro- 
cedure thus consists of fitting the previous expression with a 
pair of (Fint,T), assuming that the opacity has not significantly 
changed in the course of the observations. The antenna temper- 
ature is then given by: 
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Fig. A.l: Example of Effelsberg calibration measurements 
taken on NGC 7027 for one day at the two frequencies consid- 
ered here. Left: Flux measurements for azimuthal (black cir- 
cles) and horizontal (white circles) scans. The elevation gain 
curve is plotted, normalized to the highest flux of the graph. 
Right: fit of eq. IA.2I on data collected during a single day for 
azimuthal (full line) and horizontal (dashed lines) data. The re- 
sulting Tj and Fint are the mean of these two fits. 
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The NGC 7027 flux at 23 GHz was taken as S.llJy. We de- 
rived a conversion gain of 1 .86±0. 1 and the opacity was around 
0.2, consistent with expected atmospherical transmission in the 
centrimetric domain. Figure|A] 1 illustrates the steps involved in 
the calibration process. It shows the combined effect of antenna 
gain and opacity for data taken at various elevations. Finally, 
assuming that our sources emit in a solid angle comparable 
to the main beam, the derived fluxes are translated into an- 
tenna main beam temperatures using the antenna efficiencies 
provided by the 100-m staff (A. Kraus, priv. comm.). 
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Fig. B.l: Observed spectra (solid lines) and for the IRAM data the Gaussian fits (dashed lines) for the positions listed in Table[Tl 
Each column contains the spectra of the position labeled at the top in the following order: '^CO 1-0, '^CO 2-1, C'^O 1-0, C'^O 
2-1, C^'*S 2-1 and NH3. Missing panels were not observed. 
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Fig. B.2: Same as Figure [BTI 
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Fig. B.3: Same as Figure [BTI 



